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Abstract. Cepheid distances of local galaxies (< 7 Mpc) are 
used to study the very nearby velocity field, as pioneered by 



Sandage (1986) who also pointed out its remarkable properties: 



linearity and quietness. The new data show that the velocity 
dispersion in the distance range as seen from the barycentre of 
the Local Group 1 — 8 Mpc is as low as 38 kms^^. The local 
rate of expansion coincides with the global Hubble constant. 
Down to 1 .5 Mpc we cannot detect a deviation from the linear 
Hubble flow. This puts an upper limit for the mass of the Local 
Group, for a wide class of Friedman models, including those 
with the cosmological constant. 

Key words: Cosmology: observations - Galaxies: distances 
and redshifts - Galaxies: general - Galaxies: kinematics and 
dynamics - Local Group 



1. Introduction 



In his classical paper Sandage (1986) studied the perturbation 
of the very nearby velocity field using the available distances to 
local galaxies, mostly obtained by ground-based observations 
of brightest resolved stars and Cepheids. He concluded that the 
deviation from the linear Hubble flow possibly was detected at 
very small distances, while the quiet, linear Hubble law starts 
at about 2 Mpc. He pointed out the remarkable properties of 
the local velocity field: 1) The linearity of the velocity-distance 
relation down to small distances, 2) the closeness of the local 
and global rates of expansion, and 3) the small velocity dis- 
persion around the Hubble law. He also asked what happens 
when the distances become more accurate, in particular, will 
the velocity dispersion still decrease below cr^ — 60kms~^ 
which he derived from those data. We also note that Karachent- 



sev & Makarov ( 1996 ) found that by using both massive and 
small galaxies the velocity dispersion in the local volume is 
(T,, = 72 kms^^. 



The aim of our Letter is to use new accurate data on the dis- 
tances of nearby galaxies infeiTed for the extragalactic Cepheid 
PL-relation (cf. Sect. 2) in order to study further the behaviour 
of the Hubble law in the outskirts of the Local Group. 

2. The sample 

We make use of the the Lyon Extragalactic Cepheid Database 
(Lanoix et al. 1999b). In this database there are 32 galaxies 
with both V and I band photometry. 23 galaxies have Cepheid 
measurement from the HST and rest have groundbased mea- 
surements. The PL relation was calibrated using a sample of 
galactic Cepheids measured with the astrometric HIPPARCOS 
satellite (Lanoix et al. 1999a) and coiTected for a Malmquist- 



like bias according to Lanoix et al. (1999c) 



As we are interested in the very local velocity field we re- 
stricted our sample by requiring i?gai < 7 Mpc, where also 
the Virgo-centric flow coiTection is small. i?gai refers to the 
distance in Mpc measured from Galaxy. We were left with 14 
galaxies. The Virgo flow at larger distances was studied with 



Cepheid-based distances in Ekholm et al. (1999a) 



In Table 1 we present the data relevant for the present anal- 
ysis. In column 1 we give the PGC number and in column 2 
we give the name. In column 3 we give the distance in Mpc 
as calculated from the distance moduli provided in the Lyon 
Extragalactic Cepheid Database. In column 4 we give the dis- 
tance in Mpc as seen from the barycentre of the Local Group. In 
columns 5-7 we give the velocities: Vq is the mean heliocentric 
velocity extracted from LEDA (the Lyon-Meudon Extragalac- 
tic Database), VVahii is Vq corrected for the centroid of the 
Local Group according to Yahil (1977) and Vcon is VVahii cor- 
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rected for the virgocentric motion. We explain these corrections 
in the next section. As regards uncertainties in the parameters 
we note that the mean error for Vq is 6.4 km s~^. In Fig. |l] we 
plot the error bars reflecting the uncertainties in the Cepheid 
distances. 

Finally we remind that LMC and NGC6822 belong to the 
Milky Way subsystem and NGC598 (M33) is a member of the 
M31 subsystem. The status of the iiTegular galaxy IC1613 is 
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Table 1. In this table we present the the crucial data for the 
14 galaxies with extragalactic PL-distances. The entries are 
explained in the text. 



al. (1987 1 but difference to the coiTection of Yahil et al. (1977 ) 
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SEXB 


1.39 


1.76 


302 


139 


131 



unclear Gurzadian et al. ( |1993| ) and Rauzy & Gurza dian ( |1998| ) 
link it with the Galaxy, but van den Bergh (1999 1 recognizes 
that IC1613 may be a free-floating member of the Local Group. 
Thus it is wise to consider it as a special case in our diagrams. 

3. The method 



To study the local Hubble flow at smaller distances, it be- 
comes increasingly important to make relevant corrections to 
distances and radial velocities. The first correction, discussed 
by Sandage (1986), is due to the shift of the observer to the 
centre of expansion, which in the self-gravitating Local Group 
is not in our Galaxy, but presumably in the barycentre. By as- 
suming the mass ratio of M31 and Galaxy is Mm-h/Mq ~ 2 
Sandage set the barycentre to be on the line between M31 and 
Galaxy at 2/3 of the distance to M3 1 . This is also our choice. 
We denote the distance from the barycentre by i?2/3- 

The observed mean heliocentric velocities must also be cor- 
rected. We first correct the observed velocity to the value as it 
would be measured by an observer in our Galaxy being at rest 
relative to the centroid of the Local Group[|. This we do accord- 
ing to Yahil et al. (1977), which is also the preferred choice 



in the LEDA database. We also remind that in Ekholm et al. 



( 1999a) the correction was also made according to Richter et 



' We note that the observer at the centroid would measure a slightly 
different velocity, depending on the location of the galaxy and the na- 
ture of the velocity field close to the Local Group. For example, if the 
deviation of the velocity field from the Hubble flow is exactly spher- 
ically symmetric relative to the centroid, one may show by starting 
from the first principles of the Hubble velocity field that then the ob- 
server at the centroid would see the velocity l/yahii/ cos^. 9 is the 
angle between the centroid and the Galaxy as seen from the measured 
galaxy. In the present case, 6 is normally quite small. Hence such a 
correction makes the velocity usually only 1-5 kms^^ larger for our 
field galaxies. We prefer here to use the simpler velocity correction as 
done by Sandage (198f ). 



was quite small. The velocity correction used in reads: 
Aw = 295.4 sin Z cos b — 79.1 cos ^ cos b — 37.6 sin 6. 



(1) 



After the solar motion is thus corrected for, the infall to 
the Virgo Supercluster is the next disturbance. We correct the 
radial velocity for the virgocentric motion following Ekholm et 
al. (|l999a|): 



Vvahii ± [v{d)u - v{d)] X 



1 -sin^e/d^ + V[i^cose, 



(2) 



where (~) is valid for points closer than the tangential point 
dgai < cos and (+) for dgai > cos Q. 

O is the angular distance of a galaxy from the centre 
of Virgo and d is its distance R from the centre normal- 
ized to the distance of Virgo {d = i?/i?viigo)- Following 
Ekholm et al. ( |1999a| ) we take i?virgo ==21 Mpc and = 
57 kms~^ Mpc~^. The Hubble velocity at a distance d from 
the centre of Virgo is v((i)H = V^osm(l) x d. The cosmo- 
logical velocity of the Virgo cluster becomes V"cosm(l) = 
1200 kms^^. The predicted velocity v{d) is solved using the 



Tolman-Bondi model as described by Ekholm et al. (1999a) 



4. Results 

We present the corrected velocity vs. distance relations for the 
very local velocity field (-R2/3 < 8 Mpc) in Figs. |l] and ^ 
In Fig. |l] velocities are corrected only for the solar motion ac- 
cording to Yahil et al. (1977). We calculate the velocity dis- 



persion without the Local Group members (the first 4 sym- 
bols marked as crosses), where certainly the internal dynam- 
ics masks the possible Hubble term which is small compared 
with the value of the velocity dispersion. The velocity dis- 
persion is cr„ = 42kms^^ solved from galaxies between 
^2/3 £ [1.0,8.0] Mpc. The expected Hubble law is given as 
a solid straight line and the la dispersion as dotted lines. 

In Fig. ^ we show the influence of the correction for the vir- 
gocentric motion. Now the velocity dispersion decreases down 
to — 37 kms^^. The decrease in is not significant. This 
is expected because most of the galaxies are at large angular 
distance from the Virgo except IC4182 having <d = 26.4°. It 
is quite remarkable that after the correction this galaxy follows 
almost exactly the expected Hubble law. 

Both Figs. |l] and || support the conclusion 4 of Sandage 
( 1986| ). The observed random velocities get smaller as the un- 
certainties in the adopted distance indicator diminish. 

Fig. H shows the predictions for the velocity perturbations 
in Sandage's point-mass model for the case when the mass of 
the Local Group A/lg = 2 x IO^^Mq. For example, van den 
Bergh ( |1999|) esti mates that A/lg = 2.3 ± 0.6 x W^Mq. 
Evans et al. ( 2000| ) give an upper limit Mlg = 2.4 x lO^^M©. 

The velocities were deduced for three cosmologies using a 
numerica l algo rithm for the Tolman-Bondi calculations (Han- 
ski et al. 2000 ). The solid curve shows the prediction for the 
now preferred Friedman model Vim — 0.3 and f^A — 0.7. The 
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Fig. 1. Velocities corrected to the rest frame of the centroid of 



the Local Group according to Yahil et al. (19771. The origin 
of the distances has been changed to the centroid following 
Sandage (1986). The solid straight line is the Hubble law for 
Hq = 57 kms^^ Mpc^^. The dotted lines give the la disper- 
sion calculated for galaxies between 1 Mpc and 8 Mpc giving 
o- = 42 km Crosses refer to LMC, N6822 N224 and N598. 
The open circle is IC1613. The error bars reflect the formal 
uncertainties in the Cepheid distances. 



dashed line corresponds to the classical flat model 57^ = 1-0 
and r^A = 0.0. The dotted line is the prediction for a low mass- 
density universe — 0.1 and Ha = 0.0. flm is the mean 
density parameter of the matter and fl\ is the density parameter 
induced by the A-term. One may note that these models differ 
little from each other, which shows that the value of the ve- 
locity deflection (and the zero-velocity distance) is essentially 
determined by the mass of the LG. 

Only galaxy that shows any significant deviation from the 
expected Hubble law is IC1613 but as pointed out in Sect. 2 
its dynamical status is not well known. On the other hand M3 1 
deviates only little from the Hubble law (in fact it is within the 
la limit). 



Sandage (1986) discussed the position of the barycentre. 
He pondered whether it would be possible to reduce the scat- 
ter in the Hubble diagram by "changing the origin of the dis- 
tances to the centroid as well". We have checked the behaviour 
of the velocity dispersion a,, as a function of the position of the 
barycentre between Galaxy and M3 1 . 

Recently Evans et al. ( 200Cl| ) claimed that the halo of M3 1 is 
about as massive as the halo of Galaxy. In this case the centroid 
would be at the middle distance. Now, ay changes quite slowly 
in the region between i?i/2 and i?2/3- For example at R1/2 the 
dispersion a^ = 44 kms^^ and at -R1/4 it is ay ~ 49 kms~^. 
It is thus not possible to make any statistically relevant claims 
about the exact position of the barycentre. With a larger sample 
of field galaxies, the position of the barycentre (and hence the 
mass ratio Mm3 i/Af Galaxy) could perhaps be determined. 
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Fig. 2. As Fig. |l] except that now we have also corrected for the 
Virgocentric motions. The effect is small except for IC4182. 
This is the only galaxy that is at a small angular distance from 
Virgo cluster The velocity dispersion a^ = 37 kms~^. 
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Fig. 3. As Fig ^. Now we have added three theoretical predic- 
tions. Each of them is based on a point mass model assuming 
that the mass of the Local Group is A/lg = 2 x IO^^A/q. 
We examined three cosmologies: 1) the Friedman model with 
ri„i = 0.3 and JIa = 0.7 (the solid curve), 2) the classical flat 
model with = 1.0 and JIa = 0.0 (the dashed curve) and 3) 
a low mass density model with = 0.1 and J7a = 0.0 (the 
dotted curve). 



Finally we note that due to the lack of galaxies with 
Cepheid distances within -R2/3 = 1 Mpc and -R2/3 = 2 Mpc 
it is not possible to determine the zero-velocity surface. It is 
though important to note that galaxies within -R2/3 ~ 1 Mpc 
show practically no deflection as compared to the theoreti- 
cal curves shown in Fig. ||, except possibly the "free-floating" 
IC1613, which could fit some theoretical curve with A/lg < 



10^2 A/, 



0- 
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5. Discussion and Conclusions 
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We may conclude that all the properties of the local Hubble 
flow, mentioned by Sandage, become even more remarkable 
when more accurate distances are used. Sandage ( 1986 ) noted 
that the derived velocity dispersion around the local Hubble law 
has progressively decreased together with increasing accuracy 
in distances, from Hubble's 200 kms^^ to his 60 kms^^. Our 
results confirm that this trend continues, and ay < 40km s^^. 

The velocity-distance diagrams (Figs. ^ and |3]) show that 
the local Hubble constant is much the same as the more global 
value derived from the Tully-Fisher indicator from the KLUN 
sample (Theureau et al. 1997| ; Ekholm et al. 1999b). This also 
confirms Sandage's ( 1999| ) recent estimate that the very local 
Hubble constant is the same as the global Hq within 10 percent. 

The third important feature is the small distance where the 
Hubble law emerges. The new data show that the linear Hubble 
law extends down to at least 1.5 Mpc. Hence, within the stan- 
dard Friedman model (including the cosmological constant A; 
see Fig. one gets an upper hmit for the mass of the Local 
Group (« 2 X IO^^A/q). 

The quiet Hubble flow within the very clumpy local galaxy 
universe has always been a real riddle (Sandage et al. 1972; 



Baryshev Yu., Chemin A., Teerikorpi P., 2000, |astro-ph/0011528 
Chemin A., Teerikorpi P., B aryshev Yu., 2000, Adv, Space Res., in 



Sandage 1999). The problem of the local quiet Hubble flow was 



studied by Governato et al. (1997) using high-resolution CDM 



N-body simulations. They constructed a large sample of "Local 
Groups" and calculated the velocity dispersions in 5 Mpc vol- 
umes around the LG candidates. They found for 1 CDM 
model that the velocity dispersion is 300 — 700 km s^^ and for 
n = 0.3 CDM model 150 - 300 kms"^ They state fliat these 
simulations were unable to produce a single LG having a ve- 
locity dispersion as low as observed. 

The very local Hubble diagram offers several important ap- 
plications for cosmology when the number of accurate Cepheid 
distances to local galaxies is increased. The mass of the Local 
Group and the position of its barycenter may be determined 
m ore pre cisely. And recently it has been suggested (Chernin et 
al. ZOOCj ) that the quietness of the local Hubble flow is a sig- 
nature of the cosmological vacuum (Einstein's A-constant cor- 
responding to energy density pA — const) dominated universe 
where the velocity perturbations are adiabatically decreasing. 
For solution of both problems -the small velocity dispersion 
and the Hubble law starting i mmedi ately at 1.5 Mpc with the 
global Hq - Baryshev et al. ( 200C ) proposed that there is a 
cosmological, homogeneous dark energy (quintessence) com- 
ponent with time-variable energy density /3q(<) and equation 
of state pq — w pqc^ with w G [—1,0]. These discussions 
show that the very local volume is extremely important for the 
study of global properties of the universe. 
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